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We compute all kaon and pion parton distribution amplitudes (PDAs) to twist-three and find that 
only the pseudotensor PDA can reasonably be approximated by its conformal limit expression. 

At terrestrially accessible energy scales, the twist-two and pseudoscalar twist-three PDAs differ 
significantly from those functions commonly associated with their forms in QCD’s conformal limit. 

In all amplitudes studied, SU{3) flavour-symmetry breaking is typically a 13% effect. This scale is 
determined by nonperturbative dynamics; namely, the current-quark-mass dependence of dynamical 
chiral symmetry breaking. The heavier-quark is favoured by this distortion, for example, support is 
shifted to the s-quark in the negative kaon. It appears, therefore, that at energy scales accessible with 
existing and foreseeable facilities, one may obtain reliable expectations for experimental outcomes 
by using these “strongly dressed” PDAs in formulae for hard exclusive processes. Following this 
procedure, any discrepancies between experiment and theory will be significantly smaller than those 
produced by using the conformal-limit PDAs. Moreover, the magnitude of any disagreement will 
either be a better estimate of higher-order, higher-twist effects or provide more realistic constraints 
on the Standard Model. 

PACS numbers: 14.40.Be, 14.40.Df, 12.15.Hh, 12.38.Lg 


I. INTRODUCTION 

Knowing the behaviour of kaon and pion light-front 
parton distribution amplitudes (PDAs) is crucial in the 
analysis of a wide variety of hard exclusive processes. For 
instance, the twist-two amplitudes modulate the evolu¬ 
tion of kaon and pion elastic electromagnetic form fac¬ 
tors and both twist-two and -three amplitudes are 

critical in the study of nonleptonic R-meson decays and 
their use in constraining elements of the CKM matrix 
0. In analyses of this type, owing to the absence of re¬ 
liable computations of the PDAs, it has been common 
to employ PDAs appropriate to the conformal limit of 
QCD when estimating the size of nonperturbative fac¬ 
tors appearing in the factorisation formula for a given 
exclusive process. For example, the normalisation ex¬ 
pected for the asymptotic behaviour of the pion’s elastic 
electromagnetic form factor, FV((5^), has typically been 
based on the following conformal form for the twist-two 
PDA d,@|: = 6u(l — u). This approach to the es¬ 

timation of soft factors in hard-process factorisation for¬ 
mulae may now be reconsidered because it has become 
possible both to calculate the pointwise behaviour of me¬ 
son PDAs using continuum methods in QCD and 
to validate those results via comparisons with analyses 
©[HI of low moments of these PDAs computed using 
lattice-QCD (IQCD). 

An illustration of the impact of these advances is pro¬ 


vided by a recent study of [13 • Specifically, it has 

been shown [3 that, at all momentum scales ^ ^qcd 
which are achievable using existing or planned facilities, 
the conformal-limit twist-two PDA, provides a 

poor approximation to the light-meson twist-two ampli¬ 
tudes. Instead, the PDA is a broad, concave function, 
viz. ~ {8/TT)y/x{l — x). Evidence in suppo rt of 

this character had long been accumulating (la - fl^ but 
overlooked. Using such a dilated PDA, the soft normal¬ 
isation factor associated with the asymptotic behaviour 
of F’.n-(Q^) increases by a factor 2-3. Thus a mismatch, 
which had long appeared to be a serious discrepancy be¬ 
tween contemporary data and direct calculations on one 
hand, and the result obtained via the factorisation for¬ 
mula on the other, is transformed into near agreement. 
Hence, in what may be a significant boost to the pro¬ 
gramme [13,113: experiments at the upgraded Jefferson 
Laboratory (JLab) [13, HO] will potentially see a clear 
sign of parton model scaling for the first time in a hadron 
elastic form factor. 

It is worth highlighting that the dilation of is 

a clear expression of dynamical chiral symmetry break¬ 
ing (DCSB) on the light-front. Consequently, empiri¬ 
cal verification at JLab of the predicted normalisation of 
Ft^(Q‘^) [13 will provide novel insights into the mecha¬ 
nism that may be identified with the generation of more 
than 98% of the proton’s mass. In a related vein, it has 
been argued [3 that one of the two twist-three pion dis¬ 
tribution amplitudes, i.e. the pseudoscalar projection of 
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the pion’s light-front wave function, may be un¬ 

derstood to describe the probability distribution of the 
chiral condensate within the pion plM^. Given the 
pseudo-Goldstone-boson character of the pion and kaon, 
a comparison between the pointwise behaviour of this dis¬ 
tribution amplitude within these mesons is particularly 
interesting and potentially instructive. 


II. DISTRIBUTION AMPLITUDES AND 
BETHE-SALPETER WAVE FUNCTIONS 

A. Definitions and Observations 

A pseudoscalar meson, Pgf{q), with mass mp, pos¬ 
sesses three two-particle light-cone distribution ampli- 
tudes to twist-three, which may be expressed thus |28l | 


The behaviour of the kaon’s twist-two PDA, tpxix), is 
also determined primarily by DCSB [ij, 113 ■ Indeed, 
lPk{x) is a broad, concave and asymmetric function, 
whose peak is shifted 12-16% away from its position in 
QCD’s conformal limit. These features show that the 
heavier quark in the kaon carries more of the bound- 
state’s momentum than the lighter quark; and also that 
DCSB modulates the magnitude of flavour-symmetry 
breaking because it is markedly smaller than one might 
expect based on the difference between light-quark cur¬ 
rent masses. Combining these features of (Pk{x) with 
those of (pnix) described above, one obtains an improved 
understanding of the ratio of kaon and pion electromag¬ 
netic form factors measured at large timelike momenta 
[^ : it eliminates much of the discrepancy between ex¬ 
periment and theory which appears if the conformal-limit 
kaon and pion PDAs are used in the relevant hard scat¬ 
tering formulae [^ . 


It would thus appear that a consistent picture is emerg¬ 
ing from the confluence between continuum and lat¬ 
tice QCD studies regarding the character of light-meson 
twist-two PDAs. Namely, that at energy scales acces¬ 
sible with existing and foreseeable facilities, reliable in¬ 
sights concerning the Standard Model may only be ob¬ 
tained by using the broad, concave PDAs whose nature 
and origin we have indicated above. Herein, therefore, we 
present calculations and results for all six kaon and pion 
two-particle distribution amplitudes that appear to twist- 
three in an expansion of the light-front wave functions of 
these pseudoscalar mesons. These amplitudes have previ¬ 
ously been estimated using their properties under confor¬ 
mal transformations as the guiding principle in concert 
with QCD sum rules [2^, [2^ ; but our analysis is the first 
to use the Dyson-Schwinger equations (DSEs) [sol - l^ . 
which have both a direct connection with QCD and can 
completely chart the pointwise behaviour of these ampli¬ 
tudes. 


This document is organised as follows. In Sec.lTI] we 
introduce the distribution amplitudes and describe the 
method that will be used in their computation. Sec- 
tion|nT] provides an algebraic illustration of our tech¬ 
niques. The algebraic formulae we obtain also serve as a 
benchmark against which to evaluate the nature of our 
numerical results, described in detail in Sec. El A sum¬ 
mary and perspective are presented in Sec. El 


(0|V’/(-a;)757 • ntljg{x)\Pgf{q)) 

= fpn-q[ 

-{0\ipf{-x)ij5i^g{x)\Pgf{q)) 

pi 

= ip^P / (ib) 

{^\lpf{-x)il5(yg.vqg.ni,ipg{x)\Pgf{q)) 

= \p%n-qj\u e-“-^ («. 0 ■ (Ic) 

Here q^ = — mp0 Xg = {zj2)nfj,, with n a light-like 
four-vector, = 0, n • q = —mp; the superscript labels 
the twist-order, which will be omitted hereafter; ^ is the 
renormalisation scale; and fp, pp are, respectively, the 
pseudovector and pseudoscalar projections of the meson’s 
Bethe-Salpeter wave function onto the origin in config¬ 
uration space, explicit expressions for which are given 
in Eqs. © below. With the conventions specified by 
Eqs. ([T]), each of the PDAs is unit normalised, viz. 

du {ipp{u; C), ujp(u; (), vp(u; C)} = 1 • (2) 

The reason why we have expressed the left-hand-side 
of Eq. (fT^ in terms of a differentiated PDA will sub¬ 
sequently become apparent - see, e.g., Eq. (1^51) and 
the preceding analysis. In considering the twist-three 
amplitudes, Eqs. m and Cil), one should bear in 
mind that they are not truly independent: a three- 
particle (quark-l-antiquark-l-gluon) twist-three amplitude 
connects them [2^ . 

Note that in order to produce quantities that are gauge 
invariant for all values of x, each of the left-hand-sides in 
Eqs. (P) should also contain a Wilson line: 


yV[-x,x] = expig dafj,Ag{a), (3) 


between the quark fields. Plainly, for any light-front tra¬ 
jectory, W[—x,x\ = 1 in lightcone gauge: n • A = 0, 
and hence the Wilson line does not contribute when 


1 We use a Euclidean metric: {7/^,7,/} = 7J = 7^; 75 = 

74717273, tr[ 757 ^ 7 ,, 7 p 7 ,^] = = (i/ 2 )[ 7 p, 7 ,,]; a b = 

o-ibil and timelike => < 0. 
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this choice is employed. On the other hand, light-cone 
gauge is seldom practicable in either model calculations 
or quantitative nonperturbative analyses in continuum 
QCD. In fact, herein, as is typical in nonperturbative 
DSE studies, we employ Landau gauge because, inter 
alia it is a fixed point of the renormalisa¬ 

tion group; that gauge for which sensitivity to model- 
dependent differences between Ansdtze for the fermion- 
gauge-boson vertex are least noticeable; and a covariant 
gauge, which is readily implemented in numerical simu¬ 
lations of lattice-regularised QCD. 

We therefore proceed by assuming that W[—x,x\ is 
not quantitatively important in computation of the two- 
particle amplitudes in Eqs. O- That has been verified 
for tpp in Eq. (fTajl [s^ and it is plausible for ujp, vp. 
However, this omission should be borne in mind; and the 
validity of the assumption judged through comparisons 
with results obtained using other methods which are also 
soundly grounded in QCD. 

The value of C in Eqs. © specifies the mass-scale rel¬ 
evant to the process in which the meson is involved and 
hence at which the PDA is to be employed. The shape of 
a given PDA changes with Q hence the ^-dependence of 
the PDAs is important. The evolution equation for two- 
particle twist-two distributions is known in closed form 
and it has the solution EHi: 

oo 

ipp{u-,T) = ip^\u) ( 4 ) 

i=0,l,2,... 

= 6m(1 — u) , (5) 

where r = l/^, u = l — m, ao = l. The expansion coef¬ 
ficients {a^J"^ ,j > 1} evolve logarithmically with r: they 
vanish as r —>■ 0. These features owe to the fact that, on 
'7 ‘Aqcd — 0, QCD is invariant under the collinear con¬ 
formal group SL(2; R) [1, The Gegenbauer-a = 3/2 
polynomials correspond to irreducible representations of 
this group and hence the expansion in Eq. (|3]). 

The evolution of the two-particle twist-three distri¬ 
butions is significantly more complicated because they 
are related to the matrix elements of three-particle 
(quark-1-antiquark-|-gluon) distributions, which have non¬ 
trivial scale dependence and mix with each other under 
renormalisation [1, [H, [2^, IH, [s^. We therefore omit 
further discussion of evolution and focus instead on re¬ 
porting result computed with 

C = C 2 := 2 GeV, (6) 

a value chosen because it is typically that employed in 
numerical simulations of IQCD. 

It is worth remarking that until recently it was com¬ 
monly assumed that at any length-scale r = !/(/, an ac¬ 
curate approximation to a given PDA was obtained using 
just the first few terms in the associated conformal ex¬ 
pansion. In connection with the twist-two distribution 
in Eq. (Hall , this amounted to the assumption that an ac¬ 
curate representation of (pp{x]T) was obtained by using 


just the first few terms of the expansion in Eq. (|4|). Let us 
call this Assumption A. It has led to models for ipp{x; t) 
whose pointwise behaviour is not concave on x £ [0,^, 
e.g. to “humped” distributions [1^. Eollowing Ref. 
one can readily establish that a double-humped form for 
<Pp{x) lies within the class of distributions produced by 
a meson Bethe-Salpeter amplitude which may be charac¬ 
terised as vanishing at zero relative momentum, instead 
of peaking thereat. No ground-state pseudoscalar or vec¬ 
tor meson Bethe-Salpeter equation solution exhibits cor¬ 
responding behaviour |4l|, |43 . 

Assumption A is certainly valid on tAqcd — 0. How¬ 
ever, it is unsound at any energy scale accessible in con¬ 
temporary or foreseeable experiments. This was high¬ 
lighted in Ref. [ToI an d in Sec. 5.3 of Ref. [1^. The latter 
used the fact [4dl - [45l| that ip^^{x) can only be a good ap¬ 
proximation to a meson’s PDA when it is accurate to 
write Uv{x) « S{x), where Uv{x) is the meson’s valence- 
quark PDE, and showed that this is not valid even at en¬ 
ergy scales characteristic of the large hadron collider. An 
identical conclusion was reached in Refs. , via con¬ 

sideration of the first moment of the kaon’s PDA, which 
is a direct measure of SU (3) flavour-symmetry breaking 
and must therefore vanish in the conformal limit. Hence, 
realistic two-particle twist-two meson PDAs are necessar¬ 
ily broader than (x); in fact, much broader. It follows 
that insistence on using just a few terms in Eq. (|4]) to 
represent a hadron’s twist-two PDA will typically lead to 
unphysical oscillations, i.e. humps, just as any attempt 
to represent a box-like curve via a Fourier series will in¬ 
evitably lead to slow convergence and spurious oscilla¬ 
tions. 

The preceding observations lead one to appreciate 
that two-particle twist-three distributions evaluated at 
^2 need not closely resemble the functional forms associ¬ 
ated with their conformal limits. In our analysis of the 
PDAs in Eqs. o, we will not prejudice our analysis in 
this way. We will instead use an alternative to Assump¬ 
tion A, which is explained in the next subsection. 


B. Computing Light-front Projections of 
Bethe-Salpeter Wave Functions 

Returning now to Eqs. o, accounting for Euclidean 
metric and making use of the relationship between Bethe- 
Salpeter wave functions in configuration and momentum 
space one obtains 

fp(fp{u) = NctTZ 2 S'^{k^)j5j‘nxp{krj:kf^), (7a) 

Jdk 

/•A 

ipl,ujp{u) = NctrZ^ S'^{krj)x5Xp{krj,kf^) , (7b) 

Jdk 

1 

-Ppv'iu) = Nctr ^4 / Xp(fc^, kfj) , 

4 Jdk 

(7c) 
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where Nc = 3; the trace is over spinor indices; /j], is a 
Poincare-invariant regularisation of the four-dimensional 
integral, with A the ultraviolet regularisation mass-scale; 
^ 2 . 4 (Ci^) are, respectively, the quark wave-function and 
Lagrangian mass renormalisation constants, computed 
using a mass-independent renormalisation scheme [4^; 
and := S{n ■ — u n ■ q). 

The Bethe-Salpeter wave function in Eqs. 0 is 

Xp{kj^,kf)) = Sf{kj^)rp{kr,fj-,q)Sg{kfj), (8) 

where Pp is the Bethe-Salpeter amplitude; Sf^g are the 
dressed-quark propagators, which are usually written in 
one of the following, equivalent forms 

Sf{k) =■ kay{k'^) + ag{k'^) (9a) 

= \l[i-fkAf{k^) + Bf{k'^)] (9b) 

= Z/(fc^)/[z7 • fc-I-M/(fc^)] , (9c) 

and kr)fj = [krj -f kfj]/2, k^, = k + rjq, kg = k — {1 — ri)q, 
r] G [0,1]. Owing to Poincare covariance, no observable 
can legitimately depend on r/, i.e. the definition of the 
relative momentum. 

The pseudoscalar meson Bethe-Salpeter amplitude in 
Eq. (m has the form {£ = kgg) 


symmetry-preserving truncations of QCD’s gap and 
Bethe-Salpeter equations. That given, then, with xp in 
hand, it is straightforward to follow the procedure ex¬ 
plained in Refs. [3-[i, [2^ and thereby obtain the meson 
PDAs from Eqs. 0. The first step is to compute the 
moments 

('«a )0 = [ du{2u- l)""4>p{u) (13) 

Jo 

where (j) = ipp, wp, t>p, which are determined explicitly 
via 

= NctrZ 2 'D{n,kg,q,m)x 5 l ■ nXp{kg,kfj), (14a) 

Jdk 

[A 

= NctrZ4 ‘D{n,kg,q,m)x 5 Xpikg,kfj), (14b) 

Jdk 

\pU^a)^' 

= Nctr Z 4 / !D(n, kg, q, m) X%{kg, kg) , (14c) 

Jdk 


Tp{i-,q) = X 5 [iEp{i;q) + 7 • PFp{£;q) 

+X ■ £ Gp{£] q) + ag,i,igq,yHp{£] q)]. (10) 

Each of the scalar functions in Eq. (1101) has the following 
decomposition 


J^{£-,q)=J^o{e;q)+i-qJ^l{J;q): ( 11 ) 


with Jp,! even under {£-q) —>■ {—£-q). Herein, we consider 
the isospin-symmetric limit rriu = ^ m^. Note that 

for mesons constituted from valence-quarks with equal 
current-mass, Ti = 0 . 

At this point it is useful to expose the meaning of the 
normalisation factors in Eqs. 0 - To that end, we act 
with J? du on both sides of Eqs. (173 . (l7bl) and employ 
Eqs. 0, thereby arriving at: 


fp = NctTZ2 / ^5l‘nXp{kr):kfj), 

Jdk 

/*A 

ip^p = Nctr Z 4 / 75 Xp{kg, kg). 

Jdk 


( 12 a) 

( 12 b) 


These expressions will readily be recognised as distinct 
projections onto the origin in configuration space of 
the meson’s Bethe-Salpeter wave function, i.e. the me¬ 
son’s pseudovector and pseudoscalar decay constants 
[ 4 ^. Both express an intrinsic property of the meson 
and are equivalent order parameters for DCSB [4^. In 
this connection, the latter (pj,) has been identified as the 
in-meson chiral condensa te I2l[j25ll . 

As reviewed elsewhere |30l - l33l| , it is now possible to ob¬ 
tain realistic meson Bethe-Salpeter amplitudes by solv¬ 
ing a coupled system of integral equations; namely. 


where [n ■ qp^^‘D{n, k, q, m) = [2n ■ k-n- qp. Notably, 
beginning with an accurate form of Xp, arbitrarily many 
moments can be computed 0 
It is now useful to write 

(t>p{u) = (i)p{u) + (i)°{u), (I5a) 

</>p °(w) = {l/2)[(j>p{u) ± (j)p{u )], (15b) 

in which form the nonzero moments of (j>p{u) reproduce 
all the m-even moments of (j)p th® nonzero mo¬ 
ments of 0 p(m) are the m-odd moments of (j)p. Plainly, 
(t>p{u) = 0 for mesons comprised from valence-quarks 
with degenerate current-masses; but Eqs. m enable us 
to describe the treatment of all systems simultaneously. 

Consider now that Gegenbauer polynomials of order a, 
{C“( 2 m — 1 ) I n = 0 ,..., 00 }, are a complete orthonormal 
set on M S [ 0 , 1 ] with respect to the measure [u(l — m)]““ , 
a- = a — 1/2. They therefore enable reconstruction 
of any function defined on u G [Ojl]) and hence, with 
complete generality and to a level of accuracy defined by 
the summation upper bounds, 

(i>p'°{u)pirn 4 >p'^{.u), (16) 


^ Recall that our approach is formulated in Euclidean space. Thus, 
the moments provide a practical way to make the connection 
with Minkowski space, wherein the light-front is defined. In 
cases where the propagators and vertices take a simple form, 
direct calculation of the PDAs is straightforward. However, with 
sophisticated DSE-generated propagators and vertices, it is ad¬ 
vantageous to employ reconstruction from the moments. 
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where 

„i(j)p{u) = Na[u{l-u)\°‘^ ^ a“C“(2M-l), (17a) 

j=0,2,4,... 

^cl>${u) = N&[u{l - uT- Y. a“C'“(2u-l),(17b) 

1=1,3,... 

Na = r(2a+l)/[r(a+l/2)]^ and Cq = 1. In general, a ^ 
a because (j)p{u) and (j)p{u) are orthogonal components 
of 4>p{u). 

At this point, from a given set of 2mniax moments com¬ 
puted via Eqs. m, the even and odd component-PDAs 
may be determined independently by separately minimis¬ 
ing 

Y l(“A)m/(“A)<l - 1| > (18a) 

/ —2,4,... ,2mmax 

eZ= Y l(^k)™/(^k)0-lh (18b) 

I —1,3,... ,2?7T-n;iax 1 

over {a, 02 , 04 , ..., {a, oi, 03 , ..., where 

= [ du {2u - ly ^(t)p’^{u). (19) 

Jo 

This is the alternative to Assumption A mentioned 
above and exploited elsewhere [THIa It acknowl¬ 
edges that at all empirically accessible scales the point- 
wise profile of PDAs is determined by nonperturbative 
dynamics; and hence PDAs should be reconstructed from 
moments by using Gegenbauer polynomials of order a, 
with this order - the value of a - determined by the mo¬ 
ments themselves, not fixed beforehand. In known cases, 
involving tt-, K-, p- and ^-mesons, this procedure con¬ 
verges rapidly: jmax = 2 is sufficient [3-[l, . 


III. RESULTS: ALGEBRAIC BENCHMARKS 


In order to reliably compute moments via Eqs. (1141) . we 
follow the approach introduced in Ref. Q and develop a 
Nakanishi-like representation (sol - fs^ of the meson Bethe- 
Salpeter wave functions that are ultimately obtained via 
numerical solution of a Bethe-Salpeter equation. Before 
detailing the results of such numerical analysis, we judge 
it useful to provide an algebraic illustration of this idea. 

Consider, therefore, Eq. (im . treat f = g and intro¬ 
duce 


S{k) 


[-17 • k + M]AM(fc^), 
1 r(^ + 3/2) 
y/n r(i^-l-l) ^ 


Tp{k]q) = *75^ 




(20a) 

(20b) 

(20c) 


where Am(s) = l/[s -I- JVP] Am {s) = M^Am(s), k+^ = 
k — {1 — z)qj2. (Equations (|14a[L (jl4bp have been anal¬ 
ysed in this way elsewhere (7|, y.) Using a Feynman 


parametrisation in the resulting expression, the three de¬ 
nominators appearing in Eq. (I14cp can be combined into 
one A:-quadratic form, raised to a power that depends lin¬ 
early on A subsequent change of variables enables one 
to isolate the d'^k integration and arrive at 


—-pp f duu^v'p{u) = -ppm j duu™" ^vp{u) 

4 Jo 4 Jo 


( 21 ) 


4A,Z4 




'au {k^ + 


T{2 + v) 

T{u) 


X f dz f dadfd g^{z) a{l — a)'' ^ m AL{z,a, [3)"^ 

J-i Jo 


( 22 ) 


where I?(z, a, (3) = (1/2)(1 + a — 2af3 — z{l — a)). 

The expression within the parentheses in Eq. is 
recognisable as the value of pp that one obtains from 
Eq. (I12bl) using Eqs. (EUl) . It follows that 


duu^vp{u) = p(l -I- p) 


:J dzj dadp <;y{z) a{l — aY ^ Al{z,a,f3Y 


T{2v + 2)T{m + v + l) 
T{u + l)T{m + 2v + 2) ’ 


(23) 


which are the moments of 

Pp^(u; v) = (««)'' ■ (24) 

A QCD-like theory is obtained with p = 1, in which 
case 


vY^{u; p = 1) = 6uu = . (25) 

It is notable that similar reasoning leads to the same 
result for the twist-two PDA Q: 

(Pp^(u; p = 1) = 6uu = <p'^*(u). (26) 


There are two statements here. First, the algebraic model 
in Eqs. (l20l) produces an identical result for both the 
twist-two and pseudotensor twist-three PDAs. Second, 
both this equality and the functional form of the identi¬ 
cal results are precisely the outcomes argued to arise in 
the conformal limit of QCD [2^. 

The model defined by Eqs. ((^01) has also been used to 
evaluate an asymptotic form for wp, with the result Q 


Wp^(m;p) = 


(1-b p)r(2-H 2 p) 
2(1-b 2p)r(p)r(2-b p) 


\uuY-^ 


1-b 


(2p- l)(p-bl). ■ 
In this case, p = 1 produces 


(27) 


a;p^(u; p = 1) = 1 -b — u). (28) 
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Although the functions = 2,4,...} are ex¬ 

pected in the expansion of (u) on the domain 
tAqcd — 0, Eq. (l28l) is not the functional form antici¬ 
pated of QCD’s conformal-limit, which is instead [2^: 

= 1. (29) 

An explanation for this difference might be found in 
the conformal expansion of ujp{u) provided in Eq. (33) 
of Ref. [ 2 ^. There, the coefficient of the term 

is proportional to /a^, a leading moment of the pion’s 
three-particle (quark-|-antiquark-|-gluon) twist-three am¬ 
plitude, which vanishes logarithmically as the renormali¬ 
sation scale, is removed to infinity. Erom this perspec¬ 
tive, the term (1/2)C'2^^^^ in Eq. (E51) and the absence of 
an analogous correction in Eq. (1251) are a statement that 
Eqs. (1201) implicitly express a mixing pattern between the 
two- and three-particle twist-three amplitudes. 

We will subsequently provide comparisons between the 
forms in Eqs. (I25D -(I^ ^ and results for all the PDAs 
computed using realistic dressed-quark propagators and 
Bethe-Salpeter amplitudes. The differences will expose 
some of the impacts of nonperturbative dynamics and/or 
violations of SU{3) flavour-symmetry in QCD’s twist-two 
and -three two-particle sectors. Crucially, as we shall see, 
1 / = 1) in Eq. (1281) and ^ = 1) in Eq. (1251) 

are the natural benchmarks for all existing, realistic stud¬ 
ies of two-particle, twist-three PDAs that make reference 
to accessible energy scales. 

IV. RESULTS: NUMERICAL COMPUTATIONS 

A. Quark propagators and meson Bethe-Salpeter 
amplitudes 

We solved the s- and u- quark gap equations and the 
kaon and pion Bethe-Salpeter equations numerically, us¬ 
ing the interaction in Ref. . The infrared composition 
of this interaction is deliberately consistent with that de¬ 
termined in modern studies of QCD’s gauge sector (53 - 
[ 5 ^; and, in the ultraviolet, it preserves the one-loop 
renormalisation group behaviour of QCD so that, e.g., 
the dressed-quark mass-functions, Ms^u{p^) in Eg. (I5cl). 
are independent of the renormalisation point, which we 
choose to be C 2 - 

In completing the gap and Bethe-Salpeter kernels 
we employ two different procedures and compare their 
results: rainbow-ladder (RL) truncation, detailed in 
App. A.l of Ref. [13, which is the most widely used 
DSE computational scheme in hadron physics, whose 
strengths and weakness are canvassed elsewhere [sol - l^ : 
and the modern DCSB-improved (DB) kernels detailed 
in App. A.2 of Ref. [b^, which are the most refined ker¬ 
nels currently available [H, Ibll - f^ . Both schemes are 
symmetry-preserving; but the latter introduces essen¬ 
tially nonperturbative DCSB effects into the kernels, 
which are omitted in RL truncation and any stepwise 


improvement thereof. The DB kernel is thus the more 
realistic. 

As detailed elsewhere [b^ , this conclusion is supported 
by the agreement emerging between the “top-down” ap¬ 
proach to determining the quark-quark interaction in 
QCD, which works toward an ab initio computation of 
the interaction via direct analysis of the gauge-sector 
gap equations, and the DB kernel determined via the 
“bottom-up” scheme, which aims to infer the interac¬ 
tion by fitting data within a well-defined truncation of 
those equations in the matter sector that are relevant to 
bound-state properties. 

The gap and Bethe-Salpeter equation solutions are ob¬ 
tained as matrix tables of numbers. Computation of the 
moments in Eqs. (jl4p is cumbersome with such input, so 
we employ algebraic parametrisations of each array to 
serve as interpolations in evaluating the moments. For 
the quark propagators, we represent av,s as meromorphic 
functions with no poles on the real fc^-axis (^ , a feature 
consistent with confinement as defined through the vio¬ 
lation of reflection positivity [sTI - fs^ . Ib3 - I^ . Concerning 
the Bethe-Salpeter amplitudes, each scalar function in 
Eg. (HOD is expressed via a Nakanishi-like representation 
[SOl - l^ . with parameters fitted to that function’s first 
two (pion) or four (kaon) (. ■ q Chebyshev moments. The 
quality of the description is illustrated via the dressed- 
quark propagator in Fig.[TJ (Details of these procedures 
are presented in Appendix lAli 

Using the interpolating spectral representations, it is 
straightforward to compute arbitrarily many moments 
of the meson PDAs via Eqs. (HI- We typically em¬ 
ploy 2TO„iax = 50. The pointwise forms of the PDAs 
are then reconstructed via the “Gegenbauer-a” proce¬ 
dure described in connection with Eqs. (m-diii) above. 
Again, the procedure converges rapidly in all cases, so 
that results obtained with jmax = 2 produce < 1%. 

B. Twist-two 

The two-particle leading-twist PDAs for the pion and 
neg ative-kaoiil were computed and discussed in Refs. 0, 

I 111 . [2b[. The results obtained with the DB kernel are 

[l-ha^®C'^°'*(M-u)], (30) 
= 0.81, - 1/2, = -0.12; and 

= m</>^(M) + J (^1) 

with the functions defined in Eqs. CZl) and 

a a . 

1.42 1.14 0.074 0.07b 0.011 ’ ^ ’ 


® The PDA for the corresponding antiparticle is obtained through 
the replacement u u = I — u. 





FIG. 2. Two-particle twist-two PDAs, computed at C, 2 '- dot- 
dashed curve (dark blue) - pion, Eq. (I30II : solid curve (black) 
- kaon, Eqs. (1^ . (l32l) : dashed curve (dark green) - (/?™°‘*®*(u), 
Eq. (I33II : and dotted curve (indigo) - asymptotic distribution, 

Eq.®. 


C. Pseudoscalar, two-particle, twist-three 


FIG. 1. Functions characterising the dressed quark propaga¬ 
tor in the DB truncation. Upper panel, u/d-quark functions, 
lyg^yip^) ~ solution (open circles and squares, respectively) 
and interpolation functions (solid and long-dashed curves, re¬ 
spectively). Lower panel, s-quark functions, (Tjy(p^). Same 
legend. 

These predictions agree with the best available estimates 
from numerical simulations of IQCD [13, ITTI. 1^. 

We depict the curves of Eqs. (I30l) - (l3^ in Fig.® and 
compare them with the asymptotic two-particle distri¬ 
bution, Eq. I®) , and also with a model result Cl [13: 

^model(y) ^ (8/7r)y'M(l - u) , (33) 

which is practically indistinguishable from the DSE pre- 
diction0 Notably, the DB prediction for the second mo¬ 
ment of the pion’s twist-two PDA is in agreement with 
earlier results from IQCD Cl) [z3 confirmed by 

a more recent analysis ca , viz. 

(^A)°f=0.25 cf. 0.24 ±0.01. (34) 

These results support observations made in the Intro¬ 
duction; namely, that at accessible energy scales, two- 
particle twist-two PDAs are broad, concave functions in 
which violations of SU{3) flavour-symmetry are modu¬ 
lated by DCSB (since they are at the level of 14%). These 
issues are discussed further in Sec. lIV Cl 


^ N.B. Whereas the scale that should be associated with the model 
analyses is poorly known, our result is computed at ^ 2 - There¬ 
fore, the agreement emphasised here suggests that one is best 
advised to associate a scale of ^2 with the model results, too. 


1. Pion - oJtt 

The pion’s two-particle twist-three PDA, a;^(u), was 
computed and discussed in Ref. 0 . In addition to playing 
an important role in the study of R-meson pionic decays 
0 , ujTr{u) may be viewed as describing the light-front 
distribution of the chiral condensate within the pion. The 
chiral-limit prediction from Ref. Q is depicted in the up¬ 
per panel of Fig.® the solid curve corresponds to 

W 7 r(u; C 2 ) = Na[uu\°'^ [1 ± a 2 C^\u - u)], (35) 

with a = V — 1/2, v = 1.05, = 0.48. This is very 

close to in Eq. (l28)l . which is described by = 1, 

02 = 1 / 2 . 

It is important to be aware that only the (£;(?) 

term in Eq. (1101) provides a nonzero contribution to the 
right-hand-side of Eq. 120 when one removes the regular- 
isation scale A —>■ oo. That is because limA_>oo ■^ 4 (C) = 

0 ; and whilst the integral of the ET^{t,q) term diverges 
with A at precisely the rate required to produce a fi¬ 
nite, nonzero, A-independent result, the terms FTr{i;q), 
GTr{(;q), HT^{t,q) provide contributions to the integral 
that are finite as A —> oo and hence disappear when mul¬ 
tiplied by a renormalisation constant which vanishes in 
this limit. 

In addition, since the integral is dominated by the 
ultraviolet behaviour of the integrand, no difference in 
Bethe-Salpeter kernels at infrared momenta can have an 
impact. Owing to the following quark-level Goldberger- 
Treiman relation [13, [23 : 

Bo{P), (36) 

a pointwise expression of Goldstone’s theorem in QGD, 
the chiral-limit prediction, Eq. (1551) . is completely deter- 













FIG. 3. Pseudoscalar two-particle, twist-three PDAs, com¬ 
puted at (^ 2 . Upper panel - pion: solid curve (black), pre¬ 
diction in Eq. (13^; dot-dashed curve (dark green), QCD sum 
rules estimate ; and dotted curve (indigo), cjp^ in Eq. (I26II . 
Lower panel - kaon: solid curve (black), DB result; dot- 
dashed curve (greenh RL result; dashed curve (red), QCD 
sum rules estimate [^; and dotted curve (blue), pion predic¬ 
tion in Eq. (|3^ . 


mined by the momentum-dependence of the scalar piece 
of the self-energy associated with the dressed-quark that 
is confined within the pion. This momentum-dependence 
is the same in all DSE truncation schemes that preserve 
the one-loop renormalisation group properties of QCD, a 
fact that was confirmed in Ref. [8| by computing uJt^{x) 
in both RL and DB truncation and verifying that the re¬ 
sults are identical. These observations suggest strongly 
that the character of our result is a model-independent 
feature of strong-coupling QCD. 

Notable, too, is the agreement evident in Fig.|3] be¬ 
tween our prediction, Eq. (l35l) . and an earlier QCD sum 
rules result [2^. (A quantitatively similar distribution 
is reported in Ref. [13].) The discrepancy near the end¬ 
points of the domain of support is understandable, given 
that just low-order moments can practically be con¬ 
strained in a sum rules analysis and such moments pos¬ 
sess little sensitivity to the behaviour of uJt^ in the neigh¬ 
bourhood of the endpoints. We judge that the gener¬ 


ally good agreement with the prediction in Eq. (1351) pro¬ 
vides strong support for the model-independent nature of 
that result. This is further emphasised by the fact that 
the estimate in Ref. improves over an earlier calcu¬ 
lation [15 and, as gauged by the L^-norm, the modern 
refinement shifts the earlier estimate toward the result in 
Eq. (1551) . Significantly, all results for differ materi¬ 
ally from ojf{u) = 1. 

It is appropriate at this point to remark that numer¬ 
ical results for uip{x) have recently been obtained using 
a light-front constituent-quark model [tJ • Those results 
are in marked disagreement with all curves depicted in 
Fig. [3] they have curvature of the opposite sign on almost 
the entire domain of support. Consequently, the quark 
model results are in conflict with a model-independent 
prediction of QCD, which is deeply rooted in DCSB. It 
is probable that this defect originates in the inability 
of constituent-quark models to voraciously express chi¬ 
ral symmetry and the pattern by which it is broken in 
QCD. There are other, kindred examples, such as the 
failure of constituent-quark models [tsI] to deliver zero as 
the chiral-limit value for the leptonic decay constant of 
excited-state pseudoscalar mesons, which is also a model- 
independent corollary of DCSB [42. ItgI - Is^I . This flaw is 
apparent in any approach that does not reliably incorpo¬ 
rate and express the nature of chiral symmetry in QCD, 
e.g. Ref. [3- 

2. Kaon - ujk 

We have computed the pseudoscalar two-particle, 
twist-three parton distribution amplitude for the 
negative-kaon using both the RL and DB truncations: for 
the reasons explained above, only the Ep=if(£;g) term 
in Eq. (1101) provides a nonzero contribution. The results 
may be quoted in the form 

Wif (u; C 2 ) = m'/'f (m) + m(t)K{u) , (37) 

with the functions defined in Eqs. (CZD and 



a 

d 

of 

a“ 

ag 

RL 

0.52 

0.51 

0.49 

0.29 

0.19 

DB 

0.52 

0.52 

0.48 

0.13 

0.08 


The predictions in Eqs. (EZl), (IMl) are associated both 
with the moments listed in Table |T] and the solid- and dot- 
dashed-curves plotted in the lower panel of Fig.|51 semi- 
quantitative agreement with a QCD sum rules estimate 
is apparent. 

There are a number of important messages that may 
be read from these results. Plainly, as with the twist-two 
amplitude, the kaon distribution is skewed in favour of 
the heavier s-quark. Here, however, since this amplitude 
describes the light-front distribution of the chiral con¬ 
densate within the hadron [5 , the impact of the current- 
quark masses is most significant at the endpoints of the 
domain of support: the PDA is suppressed on a; < 1/2 
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TABLE 1. Moments (ua = 2u — 1) of the tt and _A-meson 
pseudoscalar two-particle, twist-three PDAs at <^ 2 , computed 
using Eqs. (ESI), E3 and (1381) . We also list values obtained 
with u) = Eq. (1281) . and computed using the QCD sum 
rules estimate for ujk in Ref. [^. N.B. For the kaon, is 
our most realistic result. 


(Ua) 

m = 1 

2 

3 

4 

5 

6 

u = 

0 

0.4 

0 

0.26 

0 

0.19 

iO-jY 

0 

0.39 

0 

0.24 

0 

0.18 

. ,RL 

ujk 

0.10 

0.40 

0.070 

0.25 

0.054 

0.19 

r 

UJk 

0.044 

0.40 

0.031 

0.25 

0.024 

0.19 

tax [2^ 

0.041 

0.36 

0.038 

0.22 

0.033 

0.16 


and enhanced on x > 1/2. The distortion may be quan¬ 
tified by considering a ratio, viz. 

, ^ I 1.28 RL 

f,^duu;K(u) \ 1-12 DB ’ 

where u = I — u. When the DB result for lpk(u) 
is employed, the analogous ratio evaluates to = 
1.14. It is thus evident that the magnitude of SU{?>) 
flavour-symmetry breaking in the kaon’s pseudoscalar 
two-particle, twist-three distribution is similar to that 
in the twist-two PDA. 

The magnitude of flavour symmetry-breaking exposed 
in Eq. (l3^ may also be compared with the 15% shift in 
the peak of the kaon’s valence s-quark_parton distribution 
function, s^{x), relative to (x) and the ratio of 
neutral- and charged-kaon electromagnetic form factors 
measured in e+e“ annihilation at su = 17.4 GeV^ [s^: 
\FksKl{su)\/\Fk.k+{su)\ ~ 0.12. By way of context, it 
is notable that the ratio of s-to-u current-quark masses 
is approximately 27 [88l| . whereas the ratio of nonpertur- 
batively generated Euclidean constituent-quark masses is 
typically 1.5 and the ratio of leptonic decay constants 
Ik/ f-K ~ 1-2 [8^. Both latter quantities are equivalent 
order parameters for DCSB. 

Moreover, a DSE-based computation of leptonic de¬ 
cay constant ratios yields fsJfB = 1-2 0, in ac¬ 
cord with a recent result from unquenched lattice-QCD 
Jbs/Ib = 1.22(8) [9l|, and the same DSE framework pro¬ 
duces /sif (0)//)/^(0) = 1.21 for the ratio of i? —?► K,tt 
semileptonic transition form factors at the maximum re¬ 
coil point, a value that is typical for estimates of this 
quantity: the results in Refs. |92 h^ may be summarised 
as /^k(0)//^.(0) = 1.26(5). 

It is therefore apparent that, as with ipK{u), the 
flavour-dependence of DCSB rather than explicit chi¬ 
ral symmetry breaking is measured by the skewing of 
ujk{u): SU{3) flavour-symmetry breaking is far smaller 
than one might naively have expected because DCSB im¬ 
pacts heavily on u,d- and s-quarks. 

Looking closer at the results in the lower panel of 
Fig. El one observes that the RL PDA is more skewed 
than the DB result, viz. the RL truncation allocates a 


significantly larger fraction of the in-kaon condensate to 
its valence s-quark. This feature is also highlighted by 
comparing the RL and DB results for the moments in Ta¬ 
ble U the m = 1,3,5 RL moments are noticeably larger 
than the odd moments obtained with the DB kernel. This 
is readily understood. RL-kernels ignore DCSB in the 
quark-gluon vertex. Therefore, to describe a given body 
of phenomena, they must shift all DCSB strength into 
the infrared behaviour of the dressed-quark propagator, 
whilst nevertheless maintaining perturbative behaviour 
for This requires Ms_„(p^) to be unnaturally 

large at = 0 and then drop quickly with increasing p^, 
behaviour which influences ujk{u) via the Bethe-Salpeter 
equation. In contrast, the DB-kernel builds DCSB into 
the quark-gluon vertex and its impact is therefore shared 
between more elements of a calculation. Hence smaller 
values of = 0 ) are capable of describing the same 

body of phenomena; and these dressed-masses need fall 
less rapidly in order to reach the asymptotic limits they 
share with the RL self-energies. The DB kernel therefore 
produces a more balanced expression of DCSB’s impact 
on a meson’s Bethe-Salpeter wave function and hence the 
PDAs derived therefrom provide a more realistic expres¬ 
sion of DCSB-induced skewing: DB provides the most 
realistic result [g^. 

The lower panel of Fig. El also facilitates a comparison 
between the kaon’s pseudoscalar two-particle, twist-three 
PDA and that obtained for the pion using the same ker¬ 
nel 0. Setting the asymmetry of the kaon’s distribution 
aside, the qualitative character of the distributions is the 
same: they are both maximal at the endpoints of the 
domain of support. This is a definitive signature of the 
Goldstone boson structure of these mesons, which shows 
that the chiral condensate is primarily located in com¬ 
ponents of the pseudoscalar meson wave functions that 
express correlations with large relative momenta, a fea¬ 
ture which ensures, inter alia, that light-front longitudi¬ 
nal zero modes do not play a material role in forming the 
chiral condensate [^. 


D. Pseudotensor, two-particle twist-three 

We have also computed the distribution vp{u) defined 
via Eq. using the same technique employed for the 
other PDAs. As was the case for the pseudoscalar distri¬ 
bution amplitudes, the value of the integral in Eq. (I14cl) is 
dominated by the ultraviolet behaviour of Ep{t,q)-. the 
other elements in the pseudoscalar meson Bethe-Salpeter 
amplitude, Eq. (flUl) . play no role. Once more, therefore, 
there is no model dependence in the pion computation 
and the RL and DB results are identical: 

C 2 ) = 6mm[ 1 - 0.0048 cf'^^^(2u - 1)], (40) 

which is not meaningfully distinguishable from {u) in 
Eq. (1^ . This function is depicted in Fig.E) 

Our predictions for the kaon are also simple. Their 
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FIG. 4. Pseudotensor two-particle, twist-three PDAs, com¬ 
puted at ^ 2 : dot-dashed curve (dark blue) - pion, v-n{u) in 
Eq. (I30II : solid curve (black) - kaon in DB truncation, v^^{u) 
in Eos. (I31II . (I32II : dashed curve (dark green) - kaon in RL 
truncation v^iu) in Eos. (I31II . (I32II : and long-dashed curve 
(red) - QCD sum rules result for the kaon from Ref. [^ . We 
do not plot the asymptotic form in Eq. (1251) because it is ef¬ 
fectively indistinguishable from our prediction for the pion. 


pointwise forms are well described by: 

2 

vJ:{u) = 6uu[l + '^ajcl^^^\2u- 1)], (41) 

i=l 

where T=RL, DB and 

ai 02 

RL 0.11 -0.0035 . (42) 

DB 0.049 -0.0034 

These functions are depicted in Fig. [Hand are associated 
with the moments in Table m 

In this instance, the distortion of the distributions can 
readily be measured by the shift in peak location relative 
to that of the pion: the RL amplitude peaks at tt = 0.57 
and DB at u = 0.54. One may also use defined via 
obvious analogy with Eq. (1391) : 

= 6^^ = 1.12. (43) 

Thus, in the more realistic (DB) case, the breaking of 
SU{3) flavour-symmetry in vk is 12%, as it is within the 
other kaon PDAs calculated herein. Hence, unsurpris¬ 
ingly, the effect is once again modulated by the current- 
quark mass dependence of DCSB. 

In Fig.l^ we also display a QCD sum rules estimate 
for Vk [291 . In this case, too, the result obtained via 
that route is semi-quantitatively in agreement with our 
prediction, possessing a similar level of distortion toward 
u = 1: 9% as measured by . (The sum rules result for 
the pion distribution is symmetric around u = 1/2 but 
otherwise qualitatively similar in shape and magnitude 
to the sum-rules kaon distribution, so it is not drawn.) 


TABLE II. Moments (ua = 2u — 1) of the tt and A-meson 
pseudotensor two-particle, twist-three PDAs at ^ 2 , computed 
using Eqs. (BOl), (Sg and (|42ll. We also list values obtained 
with v_p = Vp^, Eq. (1251) . and computed using the QCD sum 
rules estimate for vk in Ref. [2^ . N.B. For the kaon, is 
our most realistic result. 


{ua) 

m = 1 

2 

3 

4 

5 

6 

V = 

0 

0.2 

0 

0.086 

0 

0.048 

Vtt 

0 

0.20 

0 

0.085 

0 

0.047 


0.065 

0.20 

0.028 

0.085 

0.015 

0.047 

vT 

0.029 

0.20 

0.013 

0.085 

0.0070 

0.047 

UK [29] 

0.030 

0.20 

0.017 

0.088 

0.011 

0.049 


V. CONCLUSION 

We described calculations of the pointwise form for all 
kaon and pion dressed-quark (two-particle) parton distri¬ 
bution amplitudes (PDAs) to twist-three. These compu¬ 
tations have become possible owing to the development 
and use of novel methods for the algebraic interpolation 
of dressed-quark propagators and meson Bethe-Salpeter 
amplitudes, which are based on the notion of generalised 
perturbation theory integral representations [ 3 , 1 ^ . 

Before providing details, we list here our three pri¬ 
mary conclusions. Namely, only the pseudotensor PDA 
can reasonably be approximated by its conformal limit. 
At any realistic energy scale, the twist-two and pseu¬ 
doscalar twist-three PDAs differ markedly from the func¬ 
tions which are commonly associated with their forms 
in QCD’s conformal limit. Moreover, in all amplitudes 
studied, SU (3) flavour-symmetry breaking is typically a 
13% effect, the scale of which is determined by nonper- 
turbative dynamics, viz. the current-quark-mass depen¬ 
dence of dynamical chiral symmetry breaking (DCSB). 
The heavier-quark is always favoured by this distortion, 
e.g. support is shifted toward u = 1 in the K~. It follows 
that kaon and pion PDAs with the properties elucidated 
herein should serve as the basis for future attempts to 
access CP violation in the Standard Model. 

Turning now to specifics, the kaon and pion twist-two 
PDAs are broad, concave functions, in which the viola¬ 
tion of SU{3) flavour-symmetry is a 12-16% effect when 
measured by the appearance of asymmetry in the kaon’s 
PDA [^. IToffi^ . . All features of these PDAs are mod¬ 
ulated by DCSB and they cannot be approximated satis¬ 
factorily by ipp{u) = 6u{l — u), the form associated with 
QCD’s conformal limit, at any energy scale achievable 
with terrestrial facilities. 

The kaon and pion dressed-quark pseudoscalar twist- 
three PDAs, ujp{u), are particularly interesting. Al¬ 
gebraic analyses based on simple input [Eqs. d^ ] pro¬ 
duce asymptotic forms for the twist-two and pseudoten¬ 
sor twist-three PDAs that both coincide with the ex¬ 
pressions anticipated of QCD’s conformal limit, </?'’*(u) = 
6 u(l — u) [Eqs. (051) . (HD)]. However, kindred analysis for 
the pseudoscalar twist-three PDA produces {u) = 
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1 + {\/2)C^^'^\2u— 1) [Eq. (|28|) ]. This form locates sig¬ 
nificant strength at the endpoints of the distribution’s do¬ 
main of support, depleting the central region, in marked 
contrast to the function which is associated with this 
PDA in QCD’s conformal limit: Wp(u) = 1. Notwith¬ 
standing this, our numerical results for are 

best understood when referred to u!p^{u) as the bench¬ 
mark [Fig. [3]. The same is true of the amplitude esti¬ 
mated when conformal invariance is used as the guiding 
principle and QCD sum rules are employed to estimate 
the relevant mass-scale parameters. Thus, in this case, 
too, if a;p(u) is truly the conformal-limit result in QCD, 
then it is irrelevant to contemporary and foreseeable ex¬ 
periments. Furthermore, as with the kaon’s twist-two 
PDA, SU{3) flavour-symmetry breaking in the dressed- 
quark pseudoscalar, twist-three PDA, u!p{u), is a 12% 
effect when measured by the mass-induced distortion of 
the PDA. 

The pseudotensor dressed-quark twist-three PDAs are 
the simplest of the quantities we considered [Fig.|4]. The 
computed pion result is almost identical to the functional 
form associated with the conformal-limit, (p‘^p{u); and the 
kaon’s PDA is a modestly asymmetrised version of (f^p{u), 
with violation of S'17(3) flavour-symmetry again at the 
level of 12 %. 

It is worth reiterating some of the advantages in using 
the Dyson-Schwinger Equation (DSE) approach in stud¬ 
ies such as this. For example, the framework preserves 
the one-loop renormalisation group behaviour of QCD, so 
that current-quark masses have a direct connection with 
the parameters in QCD’s action and the dressed-quark 
mass-functions, Ms,u(p^), are independent of the renor¬ 
malisation point. Unlike other approaches to nonpertur- 
bative phenomena in continuum QCD, the renormalisa¬ 
tion point can be fixed unambiguously, as in lattice-QCD: 
it is not a parameter to be identified with some poorly 
determined “typical hadronic scale.” Moreover, one is 
not restricted to estimating a few low-order moments of 
the PDA. In working in the continuum and computing 
Bethe-Salpeter wave functions directly, the DSEs enable 
one to deliver predictions for the pointwise behaviour of 
PDAs on the full domain u G [0,1]. Importantly, those 
predictions are parameter-free and unify a meson’s PDAs 
with a diverse range of apparently distinct phenomena. 

A coherent picture has now emerged. Modern DSE 
studies predict PDAs for light-quark mesons that are typ¬ 
ically quite different from their conformal limits and these 
differences are a clean expression of DCSB on the light 
front. Notably, where a comparison is possible, the DSE 
results are consistent with those determined via contem¬ 
porary numerical simulations of lattice-regularised QCD. 
A new paradigm thus presents itself, from which it fol¬ 
lows that at energy scales accessible with existing and 
foreseeable facilities, one may arrive at reliable expec¬ 
tations for the outcome of experiments by using these 
“strongly dressed” PDAs in formulae for hard exclusive 
processes. Following this procedure, any discrepancies 
will be significantly smaller than those produced by us¬ 


ing the conformal-limit PDAs in such formulae. More¬ 
over, the magnitude of any disagreement will either be 
a better estimate of higher-order, higher-twist effects or 
provide more realistic constraints on the Standard Model. 
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Appendix A: Interpolations of propagators and 
Bethe-Salpeter amplitudes 


Here we describe the interpolations used in our eval¬ 
uation of the moments in Eq. (1141) . There are two sets 
of results to consider; viz., those obtained in RL trun¬ 
cation and those produced by the DB truncation. The 
interaction in Ref. has one parameter := Dui be¬ 
cause with rUg = constant, light-quark observables are 
independent of the value of w G [0.4,0.6] GeV. We use 
w = 0.5 GeV. 

In RL truncation, with rUg = 0.82 GeV and renor¬ 
malisation point invariant current-quark masses rhu/d = 
6.8 MeV, rhs = 162 MeV, which correspond to the fol¬ 
lowing one-loop evolved masses = 4.7 MeV, 

= 112MeV, we obtain ttItt = 0.14GeV, /tt = 
0.093 GeV and = 0.49GeV, /^ = 0.11 GeV. 

Using the DB truncation with rUg = 0.55 GeV, we 
obtain tUtt = 0.14 GeV, rriK = O.SOGeV from renor¬ 
malisation point invariant current-quark masses rhu/d = 
4.4 MeV, nis = 90 MeV, which yield = 

3.0 MeV, = 62 MeV and produce the follow¬ 

ing values of the dressed-quark mass M„(C 2 ) = 4.3 MeV, 
Afs(C 2 ) = 89 MeV, which are in fair agreement with mod¬ 
ern lattice estimates [^ . 

In interpolating the results from either truncation, the 
dressed-quark propagators are represented as [1^ 


Sfip) = 

1=1 


Z7 • p + m:J Z7 • p + m 


f* 


(Al) 


with ^nij ^ 0 Vj, so that meromorphic func¬ 

tions with no poles on the real p^-axis, a feature con¬ 
sistent with confinement [^. We find that = 2 is 






12 


TABLE A.l. Representation parameters. Eq. m - the pair 
(x, y) represents the complex number x + iy. (Dimensioned 
quantities in GeV). 


RL 

Zl 

mi 


7712 

u 

(0.38,0.71) 

(0.71,0.22) 

(0.14,0) 

(-0.78,0.75) 

s 

(0.42,0.32) 

(0.80,0.41) 

(0.12,0) 

(-1.26,0.63) 

DB 

Zl 

mi 


m2 

u 

(0.42,0.24) 

(0.44,0.19) 

(0.13,0.07) 

(-0.76,0.60) 

s 

(0.43,0.30) 

(0.55,0.22) 

(0.12,0.11) 

(-0.83,0.42) 


TABLE A.2. Representation parameters. Eqs. (IMll-irMll 
(Dimensioned quantities in GeV). 


d 




a 

A‘ 

A“ 

DB: E 1 - Ce 

0.08 

-0.70 

1.08 

3.0 

1.41 

1.0 

F 0.55 

cl/W 

0.40 

0.0 

3.0/Ai. 

1.13 

1.0 

G -0.094 

2c^ 


0.0 

1.0/[Ab]^ 

0.79 

1.0 


TABLE A.3. Representation parameters associated with 
Eqs. (1101) . (Hill . (lAlOII . (Dimensioned quantities in GeV. 
Omitted quantities are zero or unused.) 


adequate; and the interpolation parameters are listed in 
Table lA.ll (Compared with Ref. [1^ , we have made a 
minor modification of the RL s-quark parameters, which 
slightly improves the quality of the fit but otherwise has 
no observable impact.) 

The Bethe-Salpeter amplitude for a pseudoscalar me¬ 
son is given in Eg. (fTUl) . For the pion, ^i{£;q) = 0 
in Eg. (ITT]) and it is natural to choose rj = lj2. In 


case we represent the scalar functions in 
= E,F,G) by 

Eq. (Uni) 

F{£-q)=F\£-q)+F'^{£-q), 

(A2) 

F\£]q) = djr J dz<;^i^{z) 

aEAi^i£l) 




, (A3) 

E'^{£-,q) = c'^ J dz‘;,,n{z) Aai{£1) , 

(A4) 

F'^{£-,q) = 4 j"dzq.n{z)A'^pk^Al.{£l) 

(AS) 

G'^{£-q) = cl j"dz<;,niz)AlAl.J£l) , 

(A6) 


with Aa(s) = A^Aa(s), + zi ■ q, a]^ = 1 — qe, 

ttp = 1/A'p — a_F, Oq = 1/[Aq]^ — ac- H{£;q) is small, 
has little impact, and is thus neglected. The interpolation 
parameters that fit our numerical results for the Bethe- 
Salpeter amplitudes are given in Table IA.2I They were 
obtained elsewhere Q through a least-squares fit to the 
Chebyshev moments 

= - [ dx Vl - x^F{£- q)Un{x), (A7) 

^ J-i 

with n = 0, 2, where Un{x) is an order-n Chebyshev poly¬ 
nomial of the second kind. Owing to 0(4) invariance, one 
may define x = k-q/iQ, with = 1 and q = (0,0, Q, iQ). 
Only results obtained using the DB kernel for the pion 
are truly relevant herein and we therefore only list in¬ 
terpolation parameters for that case. {N.B. The overall 
multiplicative factor resulting from canonical normalisa¬ 
tion of r,n. is not included in Table IXaI ) 


RL 

Eo 

El 

Eo 

El 

Go 

Gi 

VO 

-0.71 

0.17 

1.33 

5.62 

1.0 

-0.1 

Vi 





-0.7 


V2 

1.0 

2.0 

0.0 

0.0 

0.0 

0.0 

Uo 

1.0 

0.7 

0.42 

0.21 

0.0 

0.28 

Ui 





0.25 


103f/2 

6.83 

0.36 

0.90 

0.01 

-0.01 

0.70 

no 

5 

8 

5 

8 

10 

6 

ni 
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n2 

1 

2 

1 

2 

2 

2 

A 

1.8 

2.0 

1.5 

1.6 

2.1 

1.5 

DB 

Eo 

El 

Eo 

El 

Go 

Gi 

Vo 

-0.54 

-0.1 

-0.01 

1.6 

1.5 

3.0 

ni 

-0.7 

-0.4 

-0.7 

0.8 


3.0 

V2 

1.0 

2.0 

0.0 

0.0 

0.0 

0.0 

Uo 

1.0 

0.22 

0.56 

0.11 

-0.058 

0.12 

Ui 

-2.0 

-0.5 

-0.3 

-0.65 


-1.5 

W^U2 

2.5 

0.052 

0.39 

0.001 

0.049 

-0.60 

no 

4 

8 

4 

10 

5 

8 

m 

5 

12 

6 

12 


10 

n2 

1 

2 

1 

2 

2 

2 

A 

1.35 

1.7 

1.2 

1.45 

0.8 

1.1 


Here it is worth noting another detail associated with 
the generalised spectral representations. DSE kernels 
that preserve the one-loop renormalisation group be¬ 
haviour of QCD will necessarily generate propagators 
and Bethe-Salpeter amplitudes with a nonzero anoma¬ 
lous dimension where F labels the object concerned. 
Consequently, the spectral representation must be capa¬ 
ble of describing functions of ^ = p^/Aqqj^ that exhibit 
In-^Fj^y] behaviour for ^ ^ 1. This is readily achieved by 
noting that 

where D{s) is some function. Such a factor can be mul¬ 
tiplied into any existing spectral representation in order 
to achieve the required ultraviolet behaviour. 

In connection with the dressed-quark two-particle 
twist-three distributions considered herein, it is the 
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anomalous dimension of the dressed-quark mass-function 
that must properly be represented: jp —>■ = 12/25 

in the RL and DB kernels. Owing to Eq. (IMl) . this also 
affects the pion and kaon Bethe-Salpeter amplitudes. In 
the integrands describing q) above, we therefore in¬ 
clude a multiplicative factor 



C/iA 

da Cl (a) '' 


2ni 


da(;2{a) 


{P + ai ■ q + A|)"i 

C/zAf^ 

{P + ai - q + A 0"2 


(AlO) 


{lnAl/HP+Al]y-^-. (A9) 

An analogous procedure is followed for the kaon’s 
Bethe-Salpeter amplitude. In this instance, iFi{i]q) 0 
in Eq. (ITTI) . Hence, interpolations must be found for both 
%a{i] q). The following forms are flexible enough to al¬ 
low a satisfactory representation of the numerical solu¬ 
tions to the Bethe-Salpeter equations: 


E,{i-q) 



da Co (a) 


(I/q - I/i - I/2)Af ° 
{P + ai - q + A|)"o 


where Ci/^ (a) is obtained via Eq. (I20bl) . Eor the reasons 
described above, the multiplicative factor of Eq. (IA9I) is 
included in the Eg. 1 -integrands containing C 2 (q:). Inter¬ 
polation parameters for each function iFo.iii] q) are deter¬ 
mined via a least-squares fit to that function’s n = 0, 2 
Chebyshev moments. The resulting parameter values are 
listed in Table fSTSl {N.B. The overall multiplicative fac¬ 
tor resulting from canonical normalisation of Ex is omit¬ 
ted. The function Hp is, too, because it does not have a 
noticeable effect on our results.) 
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